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I present a new method to unveil the history of cosmic accretion and the build-up of super- 
^ ■ massive black holes in the nuclei of galaxies, based on observations of the evolving radio and 

(hard) X-ray luminosity functions of active galactic nuclei. The fundamental plane of black 
hole activity discovered by Merloni, Heinz & Di Matteo (2003), which defines a universal 

■ correlation among black hole mass (M), 2-10 keV X-ray luminosity and 5 GHz radio lumi- 
nosity is used as a mass and accretion rate estimator, provided a specific functional form for 

psj . the dependency of the X-ray luminosity on the dimensionless accretion rate m is assumed. I 

J> ' adopt the local black hole mass function as derived from the velocity dispersion (cr) distribu- 

lO \ tions of nearby galaxies coupled with the AI — a relation as a boundary condition to integrate 

0^ ■ backwards in time the continuity equation for the supermassive black holes evolution, neglect- 

' ing the role of mergers in shaping up the black hole mass function. Under the most general 

assumption that, independently on M, black hole accretion proceeds in a radiatively efficient 
way above a certain rate, and in a radiatively inefficient way below, the redshift evolution of 
the black hole mass function and the black hole accretion rate function (i.e. the distribution 
of the Eddington scaled accretion rates for objects of any given mass) are calculated self- 

■ consistently. The only tunable parameters are the overall efficiency of extracting gravitational 

energy from the accreting gas, e, and the critical ratio of the X-ray to Eddington luminosity, 

Q . -^^2-iokeV,cr/iEdd = a^cr, St wMch the transition between accretion modes takes place. For 

Y\ \ fiducial values of these parameters (e = 0.1 and = 10"'^), I found that half (^ 85%) of the 

C/3 . local black hole mass density was accumulated at redshift z < 1 (z < 3), mostly in radiatively 

■ efficient episodes of accretion. The evolution of the black hole mass function between z = 

and z ^ 3 shows clear signs of an anti-hierarchical behaviour: while the majority of the most 
massive objects (A/> 10^) were already in place at z ^ 3, lower mass ones mainly grew 
at progressively lower redshift, so that the average black hole mass increases with increas- 
ing redshift. Also, the average accretion rate decreases towards lower redshift. Consequently, 
sources in the radiatively inefficient regime of accretion only begin to dominate the comoving 
accretion energy density in the universe at z < 1 (with the exact value of z depending on Xcr), 
while at the peak of the black hole accretion rate history, radiatively efficient accretion domi- 
nates by almost an order of magnitude. I will discuss the implications of these results for the 
efficiency of accretion onto SMBH, the quasars lifetimes and duty cycles, the history of AGN 
feedback in the form of mechanical energy output and, more generally, for the cosmological 
models of structure formation in the universe. 
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1 INTRODUCTION 

In the last decade, supermassive black holes (hereafter SMBH) in 
the nuclei of galaxies have been discovered at ever increasing pace, 
and are now believed to reside in most of (perhaps all) the bulges 
of present-day galaxies (Kormendy & Richstone 1995; Richstone 
et al. 1998). Moreover, there is evidence of a correlation between 
the mass of the central black holes and either the mass and luminos- 
ity (Kormendy &. Richstone 1995; Magorrian et al. 1998; McLure 
& Dunlop 2002; Marconi & Hunt 2003) or the velocity disper- 
sion (Ferrarese 8l Merritt 2000; Gebhardt et al. 2000; Merritt & 



Ferrarese 2001; Tremaine et al. 2002) of its host bulges. This has 
led to the recognition that the formation and growth of SMBH and 
of their host galaxies are related processes, and that understanding 
their evolutionary history can provide fundamental insight into the 
theories of structure formation in the universe and into the physical 
nature of AGN feedback (Silk & Rees 1998; Fabian 1999; Kauff- 
mann & Haehnelt 2000; Umemura 2001; Granato et al. 2001; Cav- 
aliere & Vittorini 2002; Wyithe & Loeb 2003; Granato et al. 2004). 

Within this framework, the most important question that needs 
to be answered is when and how the mass currently locked up in 
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SMBH was assembled. In practice, this corresponds to asldng what 
is, at any given redshift z, the number of black holes per unit co- 
moving volume per unit (base 10) logarithm of mass, i.e. the black 
hole mass function (j>M{M,z) (hereafter BHMF). In the simplest 
case of purely accretion driven evolution (i.e. assuming mergers do 
not play an important role in shaping the BHMF), the simultane- 
ous knowledge of a black hole accretion rate distribution function 
(j}m{M, z), would completely determine the evolutionary solution, 
as the two distributions must be coupled via a continuity equation 
(Small & Blandford 1992; Marconi et al. 2004, hereafter M04), to 
be solved given the appropriate boundary conditions. Moreover, if 
black holes mainly grew by accretion, a direct link should exist 
between the quasar (QSO) and active galactic nuclei (AGN) phe- 
nomena, signatures of actively accreting phases, and nearby SMBH 
tracing directly the entire accumulated mass. In fact, previous com- 
parisons between direct estimates of the local black hole mass den- 
sity, pbh.o = Pbh(2 = 0) = M4>m{z = 0) dM, and the total 
energy density liberated by powerful quasars over the cosmic his- 
tory, as originally proposed by Soltan (1982), seem to consistently 
suggest that SMBH mainly grew while they were active (Salucci et 
al. 1999; Fabian & Iwasawa 1999; Yu & Tremaine 2002, hereafter 
YT02; Elvis, Risaliti & Zamorani 2003; Cowie et al. 2003; M04). 

From the theoretical point of view, several attempts have been 
made to link the evolution of the SMBH population to analytic and 
semi-analytic models of structure formation (Efstathiou & Rees 
1988; Haenelt & Rees 1993; Haenelt, Natarayan & Rees 1998; 
Haiman & Loeb 1998; Cattaneo, Haenelt & Rees 1999; Kauff- 
mann & Haenelt 2000; Monaco, Salucci & Danese 2000; Cava- 
liere & Vittorini 2000; Volonteri, Haardt & Madau 2003; Wyithe 
& Loeb 2003; Hatziminaoglou et al. 2003; Bromley, Somerville 
& Fabian 2004; Mahamood, Devriendt & Silk 2004). A somewhat 
different approach has been followed recently by Di Matteo et al. 
(2003), who made use of large cosmological hydrodynamical sim- 
ulations to predict the evolutionary history of supermassive black 
holes growth and activity. 

The path from cosmological CDM structure formation mod- 
els (and/or simulations) to the predicted AGN evolving luminosity 
functions, however, is dotted with uncertainties regarding a number 
of important physical processes. In general, semi-analytic modelers 
need to make specific assumptions (and introduce parameters) to 
describe gas cooling in dark matter halos, the fueling of the central 
nuclear black holes, the physics of the merger events and the nature 
of the various feedback mechanisms (from star formation and/or 
from AGN). In order to constrain some of these parameters, all the 
above models, then, need to be tested against observationally deter- 
mined signatures of SMBH growth history. The usual test-benches 
are: the evolution of the AGN luminosity functions (in any spe- 
cific band), the local black holes mass function, as derived from the 
M — a relation, the slope and intercept of the M — a relation it- 
self and the spectrum and intensity of the X-ray background light, 
known to be produced by the sum of individual AGN. 

Here I would like to propose an alternative approach which is 
capable to provide a self-consistent evolutionary picture for super- 
massive black holes. Such an approach is largely based on obserx'ed 
data and only on a minimal number of theoretical assumptions (and 
parameters). These assumptions are only needed to describe the 
physics of the innermost accretion process, through which black 
holes shine, and not the more complex interplay between growing 
SMBH and their galactic environment. This is an advantage for two 
reasons: first of all, the physics of black hole accretion is reason- 
ably well understood, both theoretically and phenomenologically; 
second, the physical properties of the innermost part of an accre- 



tion flow, where the dynamics is almost completely dominated by 
the strong gravitational field of the central black holes, should not 
depend on cosmology and redshift. 

The approach followed here is similar in spirit to those of 
Small & Blandford (1992); Salucci et al. (1999); YT02; M04. More 
specifically, I will make use of the local black hole mass function 
and of the evolution of AGN luminosity functions to constrain the 
evolutionary history of SMBH. The main novelty of the present 
work is the realization that simultaneous radio and (hard) X-ray 
observations of accreting black holes can provide tight constraints 
on both the mass and the accretion rate of an active black hole, 
through the so-called "fundamental plane" relationship for active 
black holes (Merloni, Heinz & Di Matteo 2003, hereafter MHD03). 
From this, the history of SMBH growth can be followed in detail up 
to the redshift at which reliable X-ray and radio luminosity func- 
tions can be obtained. 

The structure of the paper is the following: in section|2|l will 
describe how the knowledge of the X-ray and radio luminosity 
functions of local AGN, and of the local SMBH mass function can 
be used to obtain a complete census of the local SMBH population 
and activity distribution in the form of an accretion rate distribu- 
tion function. A necessary ingredient to perform this calculation is 
the functional form that relates the observed X-ray luminosity of 
an accreting black hole to its mass and accretion rate. In section|3| 
I will describe how it is possible to use our theoretical and phe- 
nomenological knowledge of the different modes of accretion onto 
a black hole to obtain such a relation. In section |4] the method to 
calculate the redshift evolution of both black hole mass and accre- 
tion rate functions is described, and its basic assumptions clearly 
spelled out. Section |5| is then devoted to the analysis of the most 
important results of the calculations and of the main properties of 
the evolving SMBH population between z = and 2 3. A more 
general discussion of the implication of these results is presented 
in section|6| Finally, I draw my conclusions in sectionQ 

Throughout this paper, we adopt a background cosmological 
model in accordance with the Wilkinson Microwave Anisotropy 
Probe (WMAP) experiment. The model has zero spatial curvature, a 
cosmological constant, JIa ~ 0.71 a cosmological constant Ho — 
70 km s~^, dominated by cold dark matter with Qm = 0.29 and 
fii, = 0.047 (Spergel et al. 2003). 



2 THE IMPORTANCE OF ESTIMATING THE 

CONDITIONAL RADIO/X-RAY AGN LUMINOSITY 
FUNCTION 

Large optical surveys carried out in recent years (Boyle et al. 2000; 
Fan et al. 2001; Wolf et al. 2003) probe the evolution of the QSO 
luminosity function up to high redshift, and all agree in establishing 
a strong rise in their activity from the local universe up to redshift 
2 ~ 2 and a decline above 2 ~ 3. However, inferring the proper- 
ties of the entire class of accreting supermassive black holes from 
samples selected in a single waveband can be misleading, in partic- 
ular in the case of the optical band, where obscuration effects can 
be significant. Indeed, Barger et al. (2003) have clearly shown that 
optically selected broad-line AGN and QSO form only about one 
third of the X-ray background: hard X-ray selected samples, there- 
fore, provide a more direct probe of SMBH activity (see also Catta- 
neo & Bernardi 2003). Recent works by Hosokawa (2004) and M04 
have also demonstrated that the redshift evolution of the hard X-ray 
luminosity function better describes the growth history of accreting 
supermassive black holes. Discrepancies with the results obtained 
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from optically selected QSO luminosity functions for the average 
accretion efficiency and the local black hole mass density can be 
understood by taking into account the luminosity dependence of 
both obscuration and bolometric corrections in the different bands 
(see e.g. Ueda et al. 2003; Cattaneo & Bernardi 2003; Hosokawa 
2004; M04). 

Similarly, hard X-ray emission reveals fundamental properties 
of an accreting black hole: In a recent paper (MHD03) it has been 
shown that if we define the instantaneous state of activity of a black 
hole of mass M (in units of solar masses), by the radio (at 5 GHz, 
Lb) and hard (in the 2-10 keV band, Lx) X-ray luminosity of its 
compact core, and represent such an object as a point in the three- 
dimensional space (log I/R, log Lx, log M), all black holes (either 
of stellar mass or supermassive) will lie preferentially on a plane 
(the "fundamental plane" of black hole activity), described by the 
following equation: 

login = (0.60^;5;"l)logix+(0.78^S:S9)logA^+7.33^1^^(l) 

Equation Q can be inverted to relate BH masses to observed nu- 
clear radio and X-ray luminosities: 

logM ~ 1.28(log Lr- 0.60 log Lx)- 9.34 ±1.06 (2) 
= g(logLR,logLx). 

This is an entirely empirical relation, and as such is independent on 
any accretion (or jet) physical model. It shows, however, that disc 
and jet emission from active black holes of any mass, from galac- 
tic X-ray binary sources to the most powerful quasars, are physi- 
cally and observationally correlated phenomena. Moreover, as the 
fundamental plane relationship is obeyed by the intrinsic hard X- 
ray luminosities, it is basically unaffected by absorption, and there- 
fore largely independent on the validity of any specific unification 
scheme for AGN. 

One of the consequences of this relationship is that, in an ideal 
case, the conditional radio/X-ray luminosity function of active 
black holes, i.e. the number of sources per unit co-moving volume 
per unit logarithm of radio and X-ray luminosity, $c(Lr,Lx), 
could be used to reconstruct the mass function of the underlying 
black hole population. The importance of such a possibility should 
not be underestimated. For example, it could provide an alternative 
way to study the demography of the SMBH population (at any red- 
shift) to be compared with what obtained from the M — a relation 
(which the Sloan Digital Sky Survey will largely contribute to, see 
e.g. McLure & Dunlop 2004), or with any other analysis based on 
different mass estimators (see e.g. Vestergaard 2004). In alterna- 
tive, such comparisons can be used to test the redshift evolution of 
the fundamental plane relationship itself. 

Although the future goal for solving the problem at hand 
should therefore be identified with the study of large multi- 
wavelength (X-ray and radio in particular) SMBH samples, nev- 
ertheless I will argue here that it is still possible to make some 
progress with the currently available pieces of information. In fact, 
the lack of the exact knowledge of the conditional radio/X-ray AGN 
luminosity function can be (at least partially) superseded, given the 
two separate radio, (/>r(Lr, z), and X-ray, (l>x{Lx, z), luminosity 
functions at redshift z, and an independent estimate of the black 
hole mass function, <^m(M, z) at the same redshift. By taking into 
account the fundamental plane relationship, we have that the condi- 
tional luminosity function ^'c has to satisfy the following integral 
constraints: 

/•oc 

<^x(Lx)dlogLx = / *c(Lx,LR)dlogLR (3) 

"'iR,mi„ 



(^R(LR)dlogLR = / *c(Lx,LR)dlogLx (4) 
(j)M{M)dlogM ^ vE'c(Lx,LR)dlogLRdlogLx, (5) 

J J log M<g<log M + d log A/ 

where g{Ln, Lx) is defined in equation ^2)- Iri the above formulae, 
the lower end of the X-ray and radio luminosity functions, Lx.min 
and LR.min, should be chosen in such a way as to give the same 
total number of objects with mass larger than a certain minimal 
value A/niin, as computed by integrating the BHMF. 

In the following, I will make the assumption that any func- 
tion ^c(Lx,Lh,) that satisfies equations 0, J4} and Js} can be 
regarded as the true conditional radio/X-ray luminosity function. 
To begin with, I will show in the next section how to use the above 
formalism to derive informations about the state of activity of the 
local SMBH population, in the specific form of its accretion rate 
distribution function. 



2.1 The census of the local population 

Using the integral constraints provided by eqs. <3I5> . it is possible 
to deduce a conditional radio/X-ray AGN luminosity function at 
redshift zero. In order to do that, a specific choice of the observed 
luminosity (and mass) functions needs to be done. I will adopt the 
following: 

• For the hard X-ray luminosity function (HXLF) I use the re- 
cently estimated one of Ueda et al. (2003). This is arguably the 
most complete luminosity function in the 2-10 keV spectral band, 
spanning the luminosity range of 10''^'^ — 10^®'^ erg s~^, corrected 
for absorption. 

• The 5GHz radio luminosity function (RLF) of AGN is here 
derived from the lower frequency one of Willott et al. (2001), ob- 
tained from the 3CRR, 6CE and 7CRS complete samples, assuming 
for the radio spectral index a constant value qr = 0.7 to rescale 
the luminosities. Although other determinations of the local AGN 
RLF at higher frequencies are available (see e.g. Sadler et al. 2002), 
the Willott et al. one is probably the most accurate to date up to rel- 
atively high redshifts, and this will be instrumental in studying the 
growth history of the SMBH population (see section|4}. 

• The local black hole mass function is instead estimated fol- 
lowing Aller and Richstone (2002). They derived it from the local 
luminosity functions of galaxies of different morphologies given 
in Marzke et al. (1994), together with the empirical relationships 
among total luminosity, bulge luminosity and velocity dispersion. 
However, differently from Aller & Richstone (2002), we do take 
into account the effects of the scatter in the Al — a relation^ , in the 
same way as described in Yu & Lu (2004). Although the overall ef- 
fect of considering such a scatter affects the value of the total black 
hole mass density only marginally (a dispersion of 0.27 dex results 
in an increase of the local black hole mass density of a factor of 
1.2, see YT02), nevertheless, the shape of the local black hole mass 
function is strongly affected, as demonstrated by M04. 

Starting with an initial guess for ^'c, we proceed via succes- 
sive iterations, minimizing the differences between the projections 



^ I will assume thi'oughout the paper that the black hole mass-velocity dis- 
persion relation is given by log A/rh = 8.18 -I- 4.02 log((T/200km s~^), 
as discussed in Tremaine et al. (2002). 
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of the conditional luminosity function onto the X-ray and radio lu- 
minosity axes and the observed luminosity functions, until we ob- 
tain a conditional LF that simultaneously satisfies the integral con- 
straints given by eq. J3}' '^nd Js}- 

Once such an estimate of the conditional luminosity function 
is found, it is possible to derive the local distribution of the second 
fundamental physical parameter that characterizes any active black 
hole: its accretion rate in units of Eddington luminosity. 



= 0.1 



m = eMc I I/Edd 



(6) 



(e is the accretion efficiency, see section im . Such an inversion, 
however, is model dependent, as it depends on the choice of which 
spectral energy distribution should correspond to any specific cou- 
ple of fundamental parameters M and m. In practice, it corre- 
sponds to the choice of the accretion mode of a SMBH of given 
mass and 2-10 keV luminosity. This can be done by choosing a 
specific functional form Lx = Lx{M,m), as I will describe in 
more detail in the next section. 

For the moment, suffice it to say that the main results con- 
cerning the local SMBH population are summarized by the black 
solid lines in Figures|4]and|6| showing the z — 0.1 black hole mass 
and accretion rate functions. The population of local black holes 
is dominated by sources shining, in the X-ray band, below 10~^ 
of the Eddington luminosity (assuming a 10% efficiency, see be- 
low). This is in agreement with the average value found using the 
X-ray luminosity function of Seyfert 1 galaxies by Page (2001). A 
more detailed view of the local accreting SMBH population can be 
obtained by studying the mean dimensionless accretion rate as a 
function of black hole mass, defined as: 



(m(M)) 



m{M)(l>M(M)dM 



(J) 



In figureQsuch a quantity is plotted versus the SMBH mass at 
z = 0.1. It is clear that the mean accretion rate is a strong function 
of black hole mass, with small black holes accreting at a higher rate. 
Similarly, if we define the instantaneous growth rate of a SMBH of 
mass M as M/ {M[M)), we see in figure|2|that the growth time 
is very large (about one order of magnitude larger than the Hubble 
time) for the more massive holes, a result confirmed by the 5055 
optical study of 23,000 local AGN (Heckman et al. 2004). This 
implies that the very high mass SMBH must have formed at signif- 
icantly higher redshift, as we will see in detail in section|5| where 
the redshift evolution of the SMBH population will be calculated. 

Before entering into a detailed description of the results on 
the redshift evolution of the SMBH population obtained using the 
local distributions as a boundary condition, I will discuss in the 
next section the main theoretical assumptions underlying such a 
calculation. 



3 THEORETICAL ASSUMPTIONS 
3.1 Accretion efficiency 

The accretion efficiency e, appearing in equation j6} represents the 
efficiency with which gravitational energy of the matter infalling 
onto the black hole can be extracted, regardless of it being trans- 
formed into radiation or not. It is therefore an upper limit to the 
radiative efficiency, trad, and a function of the inner boundary con- 
dition of the accretion flow only. The closest the innermost sta- 
ble circular orbit (ISCO) of the accreting gas is to the event hori- 
zon, the higher the accretion efficiency. In the classical general rel- 



Log M/Mq 

Figure 1. The mean accretion rate (in units of Eddington) (m(Af)) as a 
function of black hole mass for local {z = 0.1) accreting black holes. 

z = 0.1 




Log M/Mg 

Figure 2. The instantaneous growth rate of local supermassive black holes, 
defined as M/ {M{M)) is plotted as a function of black hole mass M. 



ativistic case of test particles, the position of the ISCO depends 
on the dimensionless angular momentum of the black holes, and 
the corresponding efficiency varies between e ~ 0.06 for non- 
rotating holes and e ~ 0.42 for maximally rotating Kerr black 
holes (Novikov & Thorne 1973 (. Recently, both YT02, comparing 
the local black hole mass function with the black hole mass den- 
sity accreted during luminous QSO phases, and Elvis, Risaliti & 
Zamorani (2002), comparing the X-ray background intensity with 
the local SMBH mass density, came to the conclusion that most su- 
permassive black holes must rapidly spinning, i.e. they must have 
accreted at an average radiative efficiency higher than the canonical 
value of 10%, so that e > erad 0.1. Such a conclusion, however, 
have been questioned by M04 on the basis of a revised estimate of 
the local SMBH mass function and of the hard X-ray luminosity 
function of AGN. 

In what follows, I will assume that the accretion efficiency is a 
constant, regardless of the nature of the specific accretion mode of 
each SMBH. The radiative efficiency, instead, will be a function of 
the accretion rate, expressed though the dependence of the observ- 
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able 2-10 keV luminosity on m and of the bolometric correction, 
as described below. 



3.2 Accretion mode transition and bolometric corrections 

The relevance of the different (theoretical) accretion modes for the 
various AGN populations is still a matter of open debate. Here I will 
follow the approach of MHD03 and try to maximize the amount of 
information on the issue by comparing active black holes of differ- 
ent masses. 

Galactic (stellar mass) black holes in X-ray binaries, either 
of transient or persistent nature, are commonly observed undergo- 
ing so-called transitions, i.e. dramatic changes in their spectral and 
variability properties (I will adopt here the terminology of McClin- 
tock & Remillard 2004, which the reader is referred to for a re- 
cent review). There are at least three well defined spectral states. 
In the low/hard state the emission is dominated by a hard X-ray 
power-law with an exponential cutoff at about few 100 keV. The 
spectrum of the thermal dominant (or high/soft) state, instead, is 
dominated by a thermal component likely originated in a standard 
Shakura &. Sunyaev (1973) accretion disc, while in the steep power- 
law (or very high) state, usually associated with a source's highest 
flux level, both a thermal and a steep power-law component sub- 
stantially contribute to the spectrum. Maccarone (2003) has shown 
that the hard-to-soft transition in these systems generally occurs at 
X-ray luminosities (in the 2-10 keV band) of about Xa ~ 5 x 10""^. 
This transition is also accompanied by a "quenching" of the steady 
radio emission observed in the low/hard state. 

Given the many similarities between the high-energy spec- 
tra of galactic black holes and AGN, a similar phenomenology 
has long been searched for in active supermassive black holes. 
The reader is referred to MHD03 for a thorough discussion on 
the scale-invariant properties of black holes coupled accretion/jet 
system, and on how to constrain theoretical accretion models for 
the different classes of objects on the basis of the observed fun- 
damental plane correlation coefficients. There we showed that, for 
black holes of any mass characterized by a: = Lx/^Edd ^0.01, 
the fundamental plane relation is consistent with the most gen- 
eral theoretical relation between radio emission, mass and accre- 
tion rate expected from synchrotron emitting jets (regardless of 
their detailed geometrical and kinematical properties), provided 
that the X-ray emitting flow is radiatively inefficient (RIAF; for a 
recent review about RIAF, see Narayan 2002 or Quataert 2003). 
In this case, Lx oc m?"^ , and the radio luminosity satisfies: 
I/R oc m^'^^M^'^* = M^'^*, i.e., Lr scales with the physical 
accretion rate only. 

On the other hand, the most luminous sources, as those falling 
into the standard definition of Quasars and broad lined AGN, must 
be accreting at an higher rate, close to the Eddington one. Also their 
spectral energy distribution, usually dominated by the so-called Big 
Blue Bump (BBB: quasi thermal UV emission most likely from an 
optically thick standard accretion disc, see e.g. Malkan 1983; Laor 
1990), indicates that above a certain critical X-ray to Eddington ra- 
tio Xci, an accretion mode transition should take place to what is 
usually described as a standard, geometrically thin and optically 
thick accretion disc (Shakura & Sunyaev 1973). More recently, 
Maccarone, Gallo & Fender (2003), analysing the same AGN sam- 
ple of MHD03, have found evidence of a connection between the 
radio-quiet AGN and X-ray binaries in the thermal dominant state. 
Also the transition luminosity, a^cr, has been found to be consistent 
with the hard-to-soft transition of galactic black holes. 

The scaling of the the hard X-ray luminosity with the accretion 



rate in this high m regime, Lx oc mf, is not straightforwardly 
predicted by the standard accretion disc theory, as the origin of the 
hard X-ray emission itself is not self-consistently predicted by the 
theory (but see, for example, Merloni 2003). On the other hand, 
observational studies of the spectral energy distributions of a large 
number of QSO and AGN and comparisons of X-ray and optical 
luminosity functions of AGN can be used to put constraints on the 
value of q I Ueda et al. 2003 I. 

In light of these facts, I will adopt here the simplest possi- 
ble functional form for the Lx/L-Edd vs. m function of a broken 
power-law, bridging the low accretion rate (radiatively inefficient) 
regime and the high accretion rate one. For radiatively efficient 
sources, we should assume, by definition, that the bolometric lu- 
minosity is simply proportional to the accretion rate, Lboi oc m. 
Then, using the fitting formula for the X-ray to bolometric luminos- 
ity ratio of M04, I obtain x oc m"'^^. Summarizing, the universal 
accretion mode function adopted here has the following scalings: 



X<Xci 
X > Xc 



(8) 



The overall normalization is found by imposing continuity at Xa 
and that in the radiatively efficient regime log x/(0.76 log m) = 
— 1.5, such that AGN at the peak of the mass distribution have X- 
ray to optical ratios consistent with observations (see e.g. Vignali, 
Brandt & Schneider 2003). 



4 THE REDSHIFT EVOLUTION OF THE BLACK HOLE 
MASS AND ACCRETION RATE FUNCTIONS 

In section im i have briefly described the general method by which 
to derive the local supermassive black holes mass and accretion 
rate distribution functions given the luminosity functions of AGN 
in both radio and X-ray bands (down to a sufficiently low luminos- 
ity in each band in order to match the integrated number densities). 
From these, the redshift evolution of the SMBH population can be 
computed integrating backwards the continuity equation that de- 
scribes SMBH evolution driven by accretion only (Small & Bland- 
ford 1992; Steed & Weinberg 2004; Hosokawa 2004): 



dcj,M{M, t) d[<t)M{M,i)- {M{M,t))] 



(9) 



dt dM 
where the mean accretion rate as a function of black hole mass and 
time, (Af) can be calculated directly from the accretion rate distri- 
bution function at time t. By setting the right hand side of equa- 
tion j9j to zero, I have implicitly assumed that mergers are unim- 
portant for the black holes growth history. A thorough examination 
of the possible consequences for SMBH growth of the inclusion of 
additional merger or direct formation terms in equation ^9) can be 
found in YT02, Hosokawa (2004) or Menou & Haiman (2004). 

In practice, starting from the BHMF and the accretion rate 
function at a given redshift z, it is possible to derive the new black 
hole mass function at redshift z + dz, (j>M(M,z + dz), by just 
subtracting the mass accreted in the time interval dt = dz{dt/dz) 
calculated according to the accretion rate function of redshift z. 
This new BHMF can then be used together with the radio and X- 
ray luminosity functions at the same redshift, (^r(Lr, z + dz) and 
(pxiLx, z + dz), to obtain the new conditional luminosity function 
{Ln , Lx , z + dz), and therefore the new accretion rate function, 
and so on. Thus, the local BHMF, as determined independently 
from the galaxy velocity dispersion distribution and the M — a re- 
lation, and the local accretion rate distribution function, as derived 
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from the conditional radio/X-ray LF and a specific accretion modes 
scenario, i.e. the function a;(m), can be used together as a boundary 
condition to integrate eq. ^9} up to the redshift where the HXLF and 
the RLE of AGN can be reliably estimated. At each redshift, I dis- 
card all SMBH whose mass has decreased below A/min ~ 10^ Mq, 
so that nothing can be said, within such a scheme, about the forma- 
tion of the SMBH seeds. Furthermore, at every redshift, the values 
of Lx,min and I/R,min (see equations |3] |4j are increased to take 
into account the loss of SMBH below our threshold. In this way, at 
every redshift the number of objects above a certain minimal mass 
is always equal to the number of radio and X-ray sources above the 
corresponding limiting luminosities. 

The final result will depend on the following assumptions: 

• The X-ray and the Radio luminosity functions of AGN, if ex- 
trapolated down to low enough luminosity, do indeed describe the 
same class of objects (SMBH), and this class of objects is assumed 
to be described by the mass function (/)m . Apart from physical argu- 
ments invoking the same inner engine for radio and X-ray selected 
AGN, there are indeed strong similarities between the cosmic evo- 
lution of radio sources and of X-ray and optically selected AGN 
(Dunlop 1998; Willott et al. 2001), a further hint of the correctness 
of this assumption. 

• The fundamental plane of black hole activity is the same at 
all redshifts, and the correlation coefficients are the same as those 
observed at redshift zero. As remarked in the introduction, this is 
a natural consequence of the fact that the fundamental plane is a 
relationship among physical quantities in the innermost region of 
the coupled accretion flow-jet system. There, the strong black hole 
gravitational field dominates, and cosmological evolution is negli- 
gible. However, it should be noted here that the radio emission from 
steep spectrum AGN, whose RLF we adopt here, may be sensi- 
tive to properties of the interstellar and intergalactic medium which 
should indeed depend on redshift. Ideally, the redshift evolution of 
the AGN radio cores luminosity function should be used, as the 
fundamental plane Q is a relationship between core luminosities 
only. 

• The evolution with redshift of both <^r and (l>x is known. This 
is certainly more accurate for the hard X-ray luminosity function, 
and in the following we adopt the luminosity dependent density 
evolution (LDDE) model of the 2-10 keV luminosity function de- 
rived by Ueda et al. (2003). The analytic approximation to the 
HXLF is described in the Appendix. 

On the other hand, the high redshift evolution of the RLF is much 
less certain. The best estimate to date is probably that by Willott et 
al. (2001). This is based on three redshift surveys of flux-limited 
samples of steep spectrum sources selected at low frequencies. 
By selecting only steep spectrum sources the authors made sure 
that the effect of strongly beamed sources (which have typically 
flat spectrum) were minimized, especially for the highest redshift 
sources. However, the fundamental plane relationship was deter- 
mined by including both flat and steep spectrum sources (but ex- 
cluding beamed objects, see MHD3), and the selection was made 
at higher frequencies. Despite this potential source of uncertainty, 
I will here adopt the Willott et al. (2001) parameterization, assum- 
ing a uniform radio spectral index or = 0.7 to rescale the fluxes. 
Previous studies of redshift evolution of the RLF of flat and steep 
spectrum sources (Dunlop & Peacock 1990) have shown that al- 
though steep spectrum sources dominate number counts, the two 
populations display similar redshift evolution. The analytical ap- 
proximation to the RLF of Willott et al. (2001) used here can be 
found in the Appendix. 
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Figure 3. Redshift evolution of the comoving black hole mass density pbh 
(top panel), of the average X-ray to Eddington ratio, (ix/iEdd) (middle 
panel), and of the average SMBH mass, in units of solai' masses (bottom 
panel). The dashed horizontal line on the middle panel marks the adopted 
value of the critical accretion rate, x^, separating radiatively inefficient ac- 
cretion from radiatively efficient one: most of the SMBH growth, therefore, 
took place during episodes of radiatively efficient accretion (see also Fig.lsl. 
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Figure 4. Redshift evolution of the SMBH mass function (BHMF), from 
redshift 3.5 till redshift 0.1. Different colors and symbols correspond to 
different redshift bins. 



• The function Lx / I/Edd ~ a;(m) does not depend on the black 
hole mass and is expressible as broken power-law: at low accretion 
rates black holes accrete in a radiatively inefficient way, and x oc 
m^'^, while at high accretion rates BH are radiatively efficient, with 
the bulk of the emission being radiated in the opticalAJV bands, so 
that X (X m}' "'^ (see section|3}. 

The whole history of supermassive black hole growth can then 
be reconstructed from the evolution of X-ray and radio AGN lu- 
minosity functions with just three free parameters: the accretion 
efficiency e, the value of the critical ratio Xa that separates the ra- 
diatively inefficient and efficient regimes, and the corresponding 
critical accretion rate or, equivalently, an X-ray bolometric correc- 
tion. 
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5 RESULTS 

The family of all possible outcomes of the calculation outlined 
above is obtained by varying parameters in the relatively nar- 
row range of physically and phenomenologically realistic values: 
0.06 < e < 0.42 and 10"'' < a::cr < 10"^. However, the uncertain- 
ties on the observed luminosity functions, the high-redshift radio 
one in particular, do not allow to put tight constraints on these pa- 
rameters. More interesting, instead, is the overall trend emerging 
from this calculation for the cosmic evolution of the SMBH pop- 
ulation and of its activity level. For this reason, here I discuss in 
detail the results of a calculation performed assuming e = 0.1 and 
Xcr ~ 10"'^, leaving a discussion of the consequences of a different 
choice of parameters to section lsTI 

Figure |3| shows, in the upper panel, the redshift evolution of 
the comoving black hole mass density from z = 3.5 to 2 = 0. The 
local black hole mass density is pbh(.z = 0) = 2.6 x 10^ Mq 
Mpc~^, a value consistent with that obtained by YT02, as already 
discussed by Alter and Richstone, (2002), while the SMBH mass 
density at redshift 3 is about six times lower: pbh(2 = 3) = 
4.5 X 10" Mq Mpc~^ . Half of the total black hole mass density was 
accumulated at redshift z < 1. This is consistent with the most re- 
cent results from X-ray Background (XRB) studies. After the deep 
Chandra and XMM surveys have revealed the redshift distribution 
of the obscured sources that most contribute to the XRB (Alexan- 
der et al. 2001; Barger et al. 2002; Mainieri et al. 2002; Rosati et 
al. 2002; Hasinger 2003; Fiore et al. 2003), the newest synthesis 
models seem to suggest that indeed a substantial fraction of the 
locally measured mass density of SMBH (maybe up to 50%) was 
accumulated in (obscured)^ low-mass AGN (with M < 10* M©) at 
z < 1 (see e.g. Gandhi & Fabian 2003, or Fabian 2003 and refer- 
ences therein). This is also consistent with the general picture for 
the BHMF evolution discussed below. 

The middle panel of Figure|3|shows the evolution of the aver- 
age X-ray to Eddington rate, defined as: 



x(j)x{z)dx 



(10) 



Here <j)x {z) is the X-ray to Eddington ratio function, representing 
the number of sources per unit co-moving volume per unit loga- 
rithm of the X-ray to Eddington ratio. It is simply related to the 
accretion rate function thorough the monotonic function of equa- 
tion (SJ: 4>x{z) — 4>m{z)drn/dx. In the same panel, the dashed 
horizontal line marks the value of a^cr. The average accretion rate 
increases by almost an order of magnitude from redshift zero un- 
til z ~ 2, where the luminosity density of the hard X-ray selected 
sources peaks, and then levels off. Supermassive black holes were 
more active in the past, in the sense that their average dimension- 
less accretion rate was higher at higher redshift. A similar con- 
clusion was drawn by Small & Blandford (1992), who adopted a 
phenomenological approach close to the one followed here, and by 
Haiman & Menou (2002) and Menci et al. (2003) in the framework 
of semi-analytic models for structure formation in cold dark matter 
universes (see discussion below, §|6j. 



It is worth emphasizing that both the fundamental plane relationship 
and the HXLF used here rely on absorption corrected 2-10 keV luminosi- 
ties, and ai'e therefore unaffected by Compton thin absorption. Therefore 
throughout the paper no distinction is made (and is possible) between ob- 
scured and unobscured sources. 




Figure 5. Comparison of the evolution of comoving black hole accre- 
tion rate densities computed according to different prescriptions. The thick 
solid line shows the results from this work, calculated for e = 0.1 and 
Xcr = 10""^. The thin triple dotted-dashed line, instead, shows the results 
derived from the hard X-ray luminosity function alone, assuming radiative 
efficiency for all sources and a fixed bolometric correction (see e.g. M04). 
Thin dashed line shows the evolution of the comoving accretion rate density 
calculated from the evolution of the QSO luminosity function of Boyle et al. 
(2000) with the bolometric correction of Elvis et al. (1994). Thick dashed 
and dot-dashed lines in the same plot represent the accretion rate density 
for sources above and below the critical rate, respectively. 



The bottom panel of Fig. |3| shows the average black hole 



mass. 



{M{z)) 



J^^^^Mcl>M{z)dM 



(11) 



where at every redshift we have taken Mmin ~ 10® M©. The anti- 
hierarchical (Granato et al. 2001, 2004; M04) nature of supermas- 
sive black holes growth is summarized by this plot, showing the 
increase of the average black hole mass with increasing redshift. 

More specifically, the evolution of the shape of the black hole 
mass function, which in turn determines the evolution of (A/), is 
shown in Figure |4| As opposed to the standard picture of hierar- 
chical mass build up of dark matter halos in CDM cosmologies, 
supermassive black holes growing by accretion between z ^ 3 and 
now have a mass function which is more and more dominated by 
largest mass objects at higher and higher redshift, at least up to 
the limit where we can trust the evolution of our parametrized X- 
ray and radio luminosity functions. As it is indeed emerging from 
the study of the QSO population of the SDSS (Vestergaard 2004; 
McLure & Dunlop 2004), most of the more massive black holes 
(M > lO'') were already in place at z ~ 3. 

The history of accretion activity is instead summarized in fig- 
ure|5| where I plot the comoving black hole accretion rate (BHAR) 
density, (thick solid line). The thick dashed and dot-dashed lines 
in the same plot represent the accretion rate density for sources 
above and below the critical rate, respectively. For comparison, 
also plotted are the evolution of the comoving accretion rate den- 
sity calculated by M04 from the hard X-ray luminosity function 
alone, assuming high radiative efficiency for all sources (thin triple 
dotted-dashed line); and the corresponding quantity, calculated in- 
stead from the evolution of the QSO luminosity function of Boyle et 
al. (2000) with the bolometric correction of Elvis et al. (1994) (thin 
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Figure 6. Redshift evolution of the SMBH accretion rate function, from 
redshift 3.5 till redshift 0.1. Different colors and symbols correspond to 
different redshift bins. 




log x^^ 

Figure 7. Shaded area is the region of the parameter space (e, Xci) which 
must be exclude, as the total black hole mass density at z = is smaller 
than the total calculated mass density accreted over cosmic time. 



dashed line). The accretion rate history is therefore dominated by 
sources with high radiative efficiency, which explains why many 
authors were indeed able to explain most of the local SMBH popu- 
lation as remnants of bright AGN phases (see e.g. YT02). However, 
due to the progressive decrease of the average accretion rate with 
time, as shown in Fig. |3| the population of accreting supermas- 
sive black holes at z < 0.5 is dominated by radiatively inefficient 
sources. For this reason, the evolution of the BHAR density from 
z = to z = 1 is less rapid than what would be inferred assum- 
ing that all growing SMBH are high m objects. The consequences 
of this fact for the history of AGN feedback are discussed in sec- 
tion]^ 

Finally, in Figure |6| the evolution of the black hole accre- 
tion rate (expressed here as X-ray to Eddington ratio, Lx/L^dd) 
function is plotted for different redshift intervals, from z = 0.1 
to z — 3.5. While the number of sources accreting at low rates 
increases monotonically with decreasing redshift, the situation is 
different for rapidly accreting objects. High accretion rate sources 
(which should correspond to QSO and bright AGNs) rapidly in- 
crease in number with increasing redshift. The cut-off redshift, 
above which the number of sources declines again, is a function 
of the typical X-ray to Eddington ratio, being lower for lower ac- 
cretion rate sources. 

The combined evolution of mass and accretion rate functions 
derived here is the cause of the strong trend observed in deep X- 
ray selected samples (Cowie et al. 2003; Hasinger 2003; Fiore et 
al. 2003), where progressively lower luminosity AGN reach their 
maximal space density at progressively lower redshifts. 

5.1 Exploring the parameter space: constraints on accretion 
efficiency 

As discussed already, recent works by YT02, Elvis et al. (2002), 
M04, have shown that a direct comparison between local SMBH 
mass density and the AGN/QSO energy density integrated over cos- 
mic time can be used to constrain accretion physics and, in particu- 
lar the mean radiative efficiency of all accreting SMBH throughout 
their history. In doing so, the obvious inequality should be satisfied, 
that the local pbh,o be always larger than (or at most equal to) the 
total mass density accreted onto active black holes. If this is not 



the case, the simplest solution to the problem is to assume a higher 
accretion (and radiative) efficiency, as indeed argued by Elvis et al. 
(2002) and YT02. 

A similar line of reasoning can of course be applied to our cal- 
culations. Following the evolution of the BHMF backwards in time, 
we can search for the values of the free parameters e and Xa for 
which the SMBH mass density becomes negative in a finite time. 
For the specific x{m) and bolometric correction adopted (see §|3j, I 
found the acceptable region is bounded below by the following em- 
pirical relation: e ~ 0.2 log(a;cr) + 0.7, which is shown in Fig.0 
with the shaded area representing the excluded region of the pa- 
rameter space. The mean radiative efficiency of SMBH that these 
calculations yield is itself a function of e and Xa, as it is shown 
in figure |8| Obviously, the higher the critical X-ray to Eddington 
rate where the transition occurs, the lower the average radiative ef- 
ficiency is for any given value of e, as a larger number of objects of 
any given Lx will be in the radiatively inefficient regime of accre- 
tion at any time. 

FiguresQand|8|seem to suggest that a relatively high average 
accretion efficiency is to be preferred, implying a non zero aver- 
age spin for the SMBH population, unless Xa is lower than 10"'^. 
It should be stressed, however, that such a constraint also depends 
crucially on the value of the local black hole mass density'', and on 
the adopted bolometric correction, and that the uncertainties in the 
HXLF and, to a larger extent, in the RLE adopted here should af- 
fect the exact determination of the accepted region of the parameter 
space. 

In any case, it is worth stressing that the qualitative behaviour 
of the BHMF evolution, being driven essentially by the evolution 
of the shape of the two luminosity functions, is not modified by 
changes of e and/or Xa. In particular, the anti-hierarchical character 
of the solution found is a robust result of the approach presented 
here. 



If Pbh,o is higher, as discussed, for example, in M04, then the radiative 
efficiency of accretion is allowed to be substantially lower without violating 
any constraint 
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Figure 8. The average radiative efficiency, trad for the evolving SMBH 
population between redshift 3.5 and 0, as a function of the accretion ef- 
ficiency e. Shown are the calculations done for three values of the criti- 
cal rate: log(xcr) = —3.5 (dashed line); —3 (fiducial case, solid line) 
and —2.5 (dot-dashed line). Only values of e that are above the line 
e ~ 0.2 log(2;cr) + 0.7 (see figurelTl are shown. 

6 DISCUSSION 

In the previous section I have shown and discussed what is possible 
to deduce about the evolution of the supermassive black hole mass 
function from the coupled evolution of HXLF and RLF. The lim- 
its of this approach are a consequence of the extreme difficulty of 
obtaining high redshift luminosity functions in the radio and X-ray 
band. Little can be said then, about the sources evolution at higher 
redshifts, where there are hints of a decline of AGN activity with 
redshift. 

Theoretically, the high z evolution should be mainly driven 
by the gravitational growth of structures, with the AGN activity 
triggered by mergers, according to the standard hierarchical pic- 
ture (Efstathiou & Rees 1988; Cavaliere & Vittorini 2000; Wyithe 
& Loeb 2003; Di Matteo et al. 2003; Menci et al. 2004). In those 
phases of early structure formation and evolution, there should be 
plenty of gas available for accretion due to the frequent galaxy 
merging events which can effectively destabilize it and make it 
reach the central SMBH sphere of influence. Black hole growth 
is essentially limited by the Eddington limit (Ihe feast, according 
to Small & Blandford 1992, or the self-limited regime, Cavaliere & 
Vittorini 2000). The shape of the black hole mass function inferred 
at 2 ~ 3 (see Fig.|4}, suggesting that higher mass black holes were 
already in place at that time, requires a very rapid growth and very 
high accretion rates in the high density peaks of the cosmic den- 
sity distribution at early times, also consistent with the trend of m 
shown in Figure|3| 

On the other hand, the rapid, anti-hierarchical, evolution of 
the SMBH population between 2 = and z ^ 2.5 has always 
been difficult to incorporate into standard CDM models for struc- 
ture formation. It is interesting, in this context, to compare the main 
qualitative, anti-hierarchical, picture emerging for growing black 
holes with the evolution of the galaxy population. As for SMBH, 
the most striking and robust evolutionary trend observed for the 
galaxy population is the rapid increase of the global star formation 
rate (SFR) between 2 = 0andz~l-2 (Lilly et al. 1996; Madau 
et al. 1996). Indeed, Di Matteo et al. (2003) and M04 have already 
shown that BHAR and SFR histories have broadly similar shapes. 



with an approximately constant ratio of about few times 10""^, and 
this has also been recently measured directly by the SDSS (Heck- 
man et al. 2004). Cowie at al. (1996) have shown that such an evo- 
lution is caused by the smooth decline with redshift (from z ^ 1 
to the present) of the rest-frame K band luminosity (an indicator of 
the total stellar mass) of galaxies that undergo rapid star formation. 
This phenomenon, called 'down-sizing' is the exact analogous of 
the anti-hierarchical growth of supermassive black holes described 
here, and has been recently confirmed by Kauffmann et al. (2004), 
by studying the environmental dependence of the relation between 
star formation and stellar mass in a large number of SDSS galaxies. 

What is the physical origin of such a common behavior? The 
combined effects of the decrease in the galaxy interaction rate in 
the era of groups and clusters formation (see e.g. Cavaliere & Vit- 
torini 2000), the expansion of the universe affecting the gas cool- 
ing efficiency (Hernquist & Springel 2003; Di Matteo et al. 2003), 
the progressive depletion of the cold gas reservoirs within galaxies 
needed to power accretion and strong feedback from both stars and 
AGN (see e.g. Wyithe & Loeb 2003; Granato et al. 2004), should 
all contribute to the observed evolution at redshift below 3, which 
can be defined as that of supply-limited accretion (Cavaliere & Vit- 
torini 2000), or the famine after the feast (Small & Blandford 1992). 
Clearly, the aim of this paper is not to explore all the above theo- 
retical issues and predictions, rather to provide a detailed picture to 
compare those predictions with. Therefore, the rest of the discus- 
sion will be devoted to two specific results on the lifetimes and duty 
cycles of active black holes and on the possible implication of the 
derived SMBH growth history for the AGN feedback evolution. 

6.1 Lifetimes and duty cycles of active black holes 

In all models that try to derive the properties of the SMBH pop- 
ulation from the observed QSO evolution, a key element is repre- 
sented by the typical quasar lifetime or by the almost equivalent 
activity duty cycle (see Martini 2003 and references therein). How- 
ever, the significance of these parameters is limited to the standard 
case in which, on the basis of an observed luminosity function in 
a specific waveband, one tries to derive the distribution of either 
BH masses or accretion rates. Usually, a constant Eddington ratio 
is assumed in this case, which implies that QSO are considered as 
on-off switches. Then, the duty cycle is simply the fraction of black 
holes active at any time, and the lifetime is the integral of the duty 
cycle over the age of the universe. 

The picture discussed here is different, in that a broad distri- 
bution of Eddington rates is not only allowed, but actually calcu- 
lated for the SMBH population at every redshift. When this is the 
case, a more meaningful definition of activity lifetime is needed. 
I follow Steed & Weinberg (2004) formulation, by first defining 
the mean Eddington rate for object of mass AIq at redshift z — 
(m(A/o , z)) and then introducing the mean accretion weighted life- 
time of a SMBH with a given mass today: 

r{Mo,z)= [ {m(Mo,z'))-^dz', (12) 

The ratio of t{Mo,z) to the Salpeter time, ts = 
eMc? /L-Edd = (e/0.1)4.5 x 10*^ yrs, gives the mean number of 
e-folds of mass growth for objects with mass Mo up to redshift z. 
The ratio of t{Mq, z) to the Hubble time tuubbic{z) ~ H{z)~^, 
instead, is a measure of the activity duty cycle of SMBH. 

As I follow here a phenomenological approach based on ob- 
served luminosity functions, and no information is therefore avail- 
able on the formation and early growth of the first black holes, it 
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Figure 9. Partial mean accretion weiglited lifetimes of SMBH with inass 
today Mo, calculated for three different redshift intervals: < z < 1 
(solid line), 1 < z < 2 (dashed Hne) and 2 < z < 3 (dot-dashed line). 
The horizontal dotted line is the Salpeter time for accretion efficiency of 
10%. The accretion weighted lifetime for BH of any given mass between 
< z < 3 is the sum of the three. 



is interesting here to calculate "partial" lifetimes in a given redshift 
interval Az = {zi,Zf): 



At{Mo,Az) 



{m(Mo,z'))-£^dz'; 



(13) 



In Figure|9| I show Ar(Afo, Az) for three redshift intervals: 
0<z<l;l<z<2 and 2 < z < 3. The accretion weighted 
lifetime for BH of any given mass between < z < 3 is of course 
just the sum of the three. The anti-hierarchical nature of mass build- 
up in actively accreting AGN and QSOs is again clearly illustrated 
by this plot. In fact, the major growth episode of a SMBH must 
coincide with the period when At > tg. This happens at 2: < 1 
for Mo < 10^-®, between redshift 1 and 2 for lO'^'^ < Mo < 10* ■^ 
and at 2 < z < 3 for lO**'^ < A/q < lO* *. Supermassive black 
holes with masses larger than Mo ~ 10*'^ today, must have ex- 
perienced their major episodes of growth at redshift higher than 3. 
Black hole of lower mass today, which are also accreting at the 
higher rates in the local universe (see section IzTl drop below an 
hypothetic seed mass (here fixed at 1O*M0, but these results do 
not depend strongly on this value) and effectively "disappear" at 
higher redshift. This reflects the obvious impossibility of working 
out the initial condition of black hole growth from the local popula- 
tion evolved backwards once an object has exponentiated its mass 
just a few times. 

It also interesting to note that the objects that dominate the 
SMBH mass function today, i.e. those in the range of masses around 
W'^'^Mq, where Mo4>m{Mo, z = 0) peaks, mainly grew around 
z ~ 1, which is when most of the X-ray background light we see 
today was emitted (Hasinger 2003). 

The ratio Ar(A/o, Az)/ At[Az), where At[Az) is the time 
elapsed in the redshift interval Az, is an indication of the "par- 
tial" duty cycle of a black hole of mass Mo today in that particular 
epoch. This is shown in figure fTol for the same redshift intervals of 

fig.ia 
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Figure 10. Partial mean accretion weighted duty-cycles of SMBH with 
mass today Mo, calculated for three different redshift intervals: < 2: < 1 
(solid line), 1 < z < 2 (dashed line) and 2 < z < 3 (dot-dashed line). 



6.2 On the kinetic energy output and the history of AGN 
feedback 

The main novelty of the constraints on black hole growth history 
presented here is the opportunity to determine both mass and ac- 
cretion rate for each and every source, thanks to the fundamental 
plane relationship. As I discussed in section|3| different modes of 
accretion are expected at different m. In fact, not just the radiative 
output of an accreting black hole scales differently with m for dif- 
ferent accretion modes, but also the total kinetic energy carried by 
the jets/outflows that are responsible for radio emission. It is pos- 
sible, therefore, to derive a parallel history of the (mostly unseen) 
mechanical power output from supermassive black holes growth. 

To derive a scaling of the jet kinetic power, Wjct, with mass 
and accretion rate for sources above and below Xcr I will proceed 
in the following way. I will fist assume that the jet kinetic power at 
injection is carried by internal energy, and I assume equipartition 
between the jet magnetic field and the total pressure in the disc (see 
Heinz & Sunyaev 2003). Then Wjct oc P^^iRl oc B^A'/^ where 
Prci represents the pressure in relativistic particles at the base of 
the jet and Rs — 2GM/(? is the Schwarzschild radius. The scal- 
ing of the magnetic field, instead, can be directly inferred from the 
inclination of the fundamental plane. Equation (1 1) of MHD03 re- 
lates the observed correlation coefficients of the fundamental plane 
relation to the slope of the electron distribution in the jet, p, the 
observable radio spectral index qr, the logarithmic derivatives of 
the magnetic field intensity with respect to mass and accretion rate 
and the index q of the Lx — m relation (see section |3)- From that 
equation we have, assuming p — 2 and qr — 0: 

dlnB 

dm 
d\nB 
dM 

where we have used the results of MHD03 ^rx = 0.6 and ^rm = 
0.78. 

Thus, we obtain the expected result that, for radiatively ineffi- 
cient flows (q ~ 2.3) the total jet kinetic power is proportional to 
the physical accretion rate only Wjct oc mM oc M (see Falcke & 
Biermann 1996; Heinz & Sunyaev 2003; Fender, Gallo & Jonker 
2003). Whether such an output dominates the energy budget (ac- 



= 0.35 g Crx 0.21 g 

= 0.35(eRX +Crm) - 1 -0.51, 



(14) 
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Figure 11. Redshift evolution of the total integrated mechanical power 
Wjct/c^ from accreting black holes per unit co-moving Mpc (triple dotted- 
dashed line). For comparison, the evolution of the total black hole accretion 
rate density (solid line) is shown together with the two separate contribu- 
tions from the sources accreting above(dashed line) and below (dot-dashed 
line) the critical rate, all taken from Fig.lsl 



cording to the so-called ADIOS picture, Blandford & Begelman 
1999, and as proposed also by Fender, Gallo & Jonker 2003) or not 
depends however on the dynamics of the innermost disc-jet cou- 
pling and, from the observationally point of view, on the radiative 
efficiency of the jets. On the other hand, SMBH accreting above 
the critical rate x^-c, for which we have assumed q — 0.76, should 
obey the scaling Wjct oc Mm°-^^. 

Given the above scaling relations for Wjct, it is possible to 
calculate the mechanical power output from each accreting black 
hole and the total integrated one, as a function of redshift. For the 
sake of simplicity I assume that indeed the total kinetic power of the 
jet/outflow from radiatively inefficient black holes dominates as a 
sink of energy. Therefore, the calculated Wjet should be considered 
an absolute upper limit to the jet/outflow kinetic power. The results 
are shown in Figure fTTI 

Because the relative contribution of the mechanical energy 
output is much larger in low accretion rate sources, the history of 
Wjet does not follow that of the most luminous sources, but instead 
exhibit a much weaker evolution, both at low and high redshifts. 
Also, because the number of low accretion rate sources increase 
with time (due to the overall decrease of the average accretion rate, 
see previous section), the peak of the mechanical power output lies 
at a much lower redshift than the peak of the radiative energy out- 
put history of AGN, which should have interesting implications for 
the formation and dynamical evolution of clusters of galaxies. 



7 CONCLUSIONS 



ous attempts to study the evolution of the supermassive black holes 
population by looking at the evolution of a single AGN luminosity 
function. 

Here, the redshift evolution of the SMBH mass function be- 
tween z — and z ~ 3 is evaluated integrating backwards in time 
a continuity equation for the black hole population in which the role 
of mergers is neglected. The local black hole mass function, esti- 
mated by applying the correlation between black hole mass and ve- 
locity dispersion to the velocity distribution function of local galax- 
ies, is used as a boundary condition for the continuity equation. The 
solution to this equation is uniquely determined once the accretion 
efficiency is specified, together with the function Lx/ LEddirn), 
that links the observed Eddington scaled hard X-ray luminosity of 
an accreting black hole to its accretion rate. For the latter, I have 
assumed a very general form of a double power-law, correspond- 
ing to a radiatively inefficient regime at low accretion rates, and a 
radiatively efficient one at high m. 

By comparing the local observed BHMF with the total mass 
accreted onto SMBH over their history, it is possible to put simul- 
taneous constraints on the accretion efficiency and on the critical 
value of the accretion rate, Xct, at which the transition takes place 
between the two accretion modes. 

The main results of this work are the following. For fiducial 
values of the parameters (e = 0.1 and Xa = 10""^), half (~ 85%) 
of the local black hole mass density was accumulated at redshift 
z < 1 (z < 3), mostly in radiatively efficient episodes of accretion. 
Qualitatively (i.e. independently on the values of these two param- 
eters), the evolution of the black hole mass function between z = 
and 2 ~ 3 shows clear signs of an anti-hierarchical behaviour. This 
is a purely phenomenological assessment, and reflects the fact that, 
while the majority of the most massive objects (M > 10^) were 
already in place at z; ~ 3, lower mass ones mainly grew at pro- 
gressively lower redshift, so that the average black hole mass in- 
creases with increasing redshift. On the other hand, the average ac- 
cretion rate decreases towards lower redshift. Therefore, sources in 
the RIAF regime of accretion only begin to dominate the comoving 
accretion energy density in the universe at jz < 1 (with the exact 
value of z depending on a::cr), while at the peak of the black hole 
accretion rate history, radiatively efficient accretion dominates by 
almost an order of magnitude. By carefully evaluating the contribu- 
tions to the total black hole accretion energy density from the differ- 
ent modes of accretion as a function of redshift, I show how to de- 
rive a physically motivated AGN feedback history. This, as well as 
the anti-hierarchical behaviour of SMBH growth described above, 
may be of some importance for cosmological models of structure 
formation in the universe. 
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I have presented a new method to study the growth of accreting su- 
permassive black holes, based on the simultaneous evolution of the 
AGN radio and hard (2-10 keV) X-ray luminosity functions. The 
method is based on the locally observed trivariate correlation be- 
tween black hole mass. X-ray and radio luminosity (the so-called 
fundamental plane of black hole activity, MHD03). Thanks to this 
correlation, it is possible for the first time to break the degeneracy 
between luminosity, mass and accretion rate that affected all previ- 
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APPENDIX A: ANALYTICAL APPROXIMATION OF THE 
X-RAY LUMINOSITY FUNCTION 

I adopt here the functional form for the HXLF described in Ueda et 
al. (2003), which has the following analytic approximation: 



(Al) 



where the local X-ray luminosity function is expressed as smoothly 
connected double power-law: 



(A2) 



while the evolutionary part is expressed as 

n(z L)-^ + z<z,{Lx) 
''^^'^''> \ r!{z,)[{l + z)/{l + z,{Lx))r z>z,{Lx) ^ ' 

and 



The anti-hierarchical growth of supermassive black holes 

For a ACDM universe, the best fit parameters are: A = 5.04 ± 

0.33 X 10-'^/i?o Mpc-^ L, = lo(*^'®''-" '^)ftTO ergs-^7i = 
0.86 ± 0.15; 72 = 2.23 ± 0.13; pi = 4.23 ± 0.39; p2 = -1.5 
(fixed); z* = 1.9 (fixed); La = lO^^'^/ifo^ erg s"^ (fixed); a = 
0.335 ± 0.070. 



APPENDIX B: ANALYTICAL APPROXIMATION OF THE 
RADIO LUMINOSITY FUNCTION 

The 5GHz RLF is derived from the low-frequency (151 MHz) one 
of Willott et al. (2001), assuming a fixed radio spectral index of 
qr = 0.7 to rescale the luminosities. 

The analytical approximation is that of model C of Willott et 
al.(2001), and is given by a sum of two differently evolving popu- 
lations (a low and a high luminosity one): 



t>R,15l{Ll51,z) = <pn,\ +^^R,h- 
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These two populations evolve differently with redshift, ac- 
cording to the following expressions: 
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The high luminosity redshift evolution fhiz) is given by: 
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The best-fitting parameters to the observed luminosity func- 
tion are: \og4>'l^i = -7.121°;;;;, ai = 0.539 ± 0.02, logLi,. = 



26.10_,, („,, 210 



-6.20 



+0.09 



ah 



1.91 ± 0.16, zu 



-0.11' oi 

0.71 ± 0.10, h = 
2.271°;;?, logLh,, 
= 0.56 ± 0.05, 2h2 = 



4.3oi;;;^I, iog<AR.h = 
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